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The frequency uj and temperature T dependent complex conductivity a of two weakly coupled 2D XY models 
subject to the RSJ dynamics is studied through computer simulations. A double dissipation-peak structure in 
Re[a><x] is found as a function of T for a fixed frequency. The characteristics of this double-peak structure, as 
well as its frequency dependence, is investigated with respect to the difference in the critical temperatures of the 
two XY models, originating from their different coupling strengths. The similarity with the experimental data in 
Festin et al. [Physica C 369, 295 (2002)] for a thin YBCO film is pointed out and some possible implications are 
suggested. 



The thermally excited two-dimensional (2D) 
vortex fluctuations drive the phase transition be- 
tween the superconducting and normal state in 
many 2D systems. Such transition is of the 
Kosterlitz-Thouless (KT) type: below the KT 
transition temperature T c vortices are bound in 
pairs with total vorticity zero, and as the tem- 
perature is increased across T c from below these 
vortex pairs start to unbind |5J]. It means that 
the dominant characteristic physical features in a 
region close to the KT transition are associated 
with vortex pair fluctuations. One of the physical 
quantities which contain the information about 
the feature of vortex dynamics is the frequency uj 
dependent complex conductivity <j(uj) of the sam- 
ple |,|. In the presence of 2D fluctuation effects 
a(uj) can be expressed as <t(uj) = — p {T) / iuje(u) , 
where l/e(uj) is the dynamic dielectric function 
which describes the effect of pair motion, and 
Po(T) is the bare superfluid density J^fl. The 
measurements of the superconducting transition 
by means of cr(u>) in a typical experiment for 
a fixed frequency should show a single peak in 
Rc[cjct(o;)] at a frequency dependent temperature 
T u . This peak represents dissipation losses while 
a sharp decrease in Im[-wcr(w)] at the same tem- 
perature is a consequence of the loss of a super- 
fluid response. Indeed, such behaviors of <j{uj) 
have been confirmed in many experiments on 2D 
superconductors, as well as on high-T c supercon- 
ductors |,|. 



The present investigation is inspired by the re- 
cent experimental results obtained by Festin et 
al. U for a 1500 A thin YBCO film: A very strik- 
ing double peak structure in Re[cjcr(w)] is found; 
two rapid drops of Im[— ujcf(uj)] at different T are 
observed. The data by Festin et al. are repro- 
duced in Figs. |l| and @. One possible explanation 
for the double peak is that the epitaxial grown 
YBCO film is split into an upper and lower part 
due to a slightly different oxygen contents 0. 
From this perspective the sample would consist 
of two weakly coupled parallel superconducting 
parts. 

In order to investigate this scenario we use the 
two weakly coupled 2D XY models, the Hamil- 
tonian of which is written as 
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where Hi and H2 are the usual 2D XY Hamilto- 
nians with the coupling strengths Ji and J2 for 
the first (lower) and the second (upper) planes, 
respectively, the summation Y]/Aj\ is over all near- 
est neighbor pairs in each plane, and H mt with 
the coupling strength Ji describes the coupling 
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between the planes. To study the dynamics of 
the system, we use the equations of motion of 
the standard resistively-shunted junction (RSJ) 
dynamics subject to the periodic boundary con- 
dition and integrate the equations of motion us- 
ing the second-order algorithm with the time step 
At = 0.05. We also apply the fast Fourier trans- 
formation method to speed up the calculations 
(see e.g., Ref. || for details). 

The 2D XY model on the square lattice under- 
goes a KT transition at T c w 0.89 J H, where J 
is the Josephson coupling strength. Accordingly, 
two planes with different coupling constants then 
undergo two separate phase transitions at the dif- 
ferent temperature T c w 0.89Ji and 0.89J2 when 
the interplane coupling vanishes. The main out- 
put from the simulation are the dynamic dielec- 
tric function l/e(w), and the helicity modulus 7, 
which basically measures the stiffness of the sys- 
tem to a twist in the phase of the order parameter 
and is proportional to the renormalized superfluid 
density p — p o /e(0) Qj. From a knowledge of 
l/e(u>) and p it is straightforward to analyze the 
behavior of the conductivity a(u>). 

In order to get realistic parameters for the two 
weakly coupled 2D XY models in connection with 
the experimental result by Festin et al. one would 
like to have some reasonable estimates of the cou- 
pling constants J\, J2, and J±. We estimate J\ 
and J2 from the experimental data reproduced 
in Fig. [|. In order to do it we use the rela- 
tion lim^-jo Im[— u>a(tu)] = po/e(Q) and note that 
1/e r* 1 just below the transition and is just 
above. Thus po oc J jlj] may be roughly estimated 
by the two heights of the rapid drops in Fig. |l| for 
the curve corresponding to the smallest frequency. 
From this we get the ratio J2/J1 ~ 0.5 (ratio be- 
tween the dashed vertical lines in Fig. Q). The 
corresponding curves for Re[cJ<7(w)] are shown 
in Fig. ||, where the characteristic double peaks 
at different temperatures are clearly exhibited. 
J± can be estimated from the knowledge of the 
anisotropy parameter T defined as T = \f J\/ J±_ 
and found to be equal to 7 for YBCO 

First we present the result for our model with- 
out coupling between the two planes. Jj_ = 
corresponds to the case of no supercurrent flow- 
ing between the planes. Figures || and || show the 




89.4 89.6 89.8 90 90.2 90.4 



T 

Figure 1. Experimentally mea- 

sured Im[— uhj(lo)] vs T for uj = 
17 mHz, 170 mHz, 1.7 Hz, 17 Hz, 170 Hz, 1.7 kHz, 
170 kHz (from the left to the right). Horizontal 
lines shows respective zero levels for transi- 
tions and vertical lines indicates the heights of 
the rapid drop of Im[-wff(w)] which gives a 
rough measure of the superconducting electrons 
involved. The data are taken from Ref. ||. 



behavior of Re[wcr(w)] and Im[— ujcf(uj)] obtained 
from the simulations with J\ = 1, J2 = 0.5 and 
J± = 0. The similarities to the experimental re- 
sults in Figs, [l] and || are striking: Re[wcr(w)] in 
Fig. |] again displays two distinct peaks at differ- 
ent temperatures while Im[— u>a(uj)\ in Fig. |^ has 
two regions with increased drops as a function of 
T. Dashed line in Fig. |^ represents the behavior 
of 7 or equally the behavior of Im[— wcr(w)] in the 
limit uj — » which shows one drop at T c \ ea 0.89 
and another at T C 2 « 0.45. These similarities 
between experiments and our simulations for the 
coupled 2D XY models are further substantiated 
when one compares the frequency dependences: 
As the frequency is increased both T c \ and T C 2 
increase as is reflected in the positions of peaks 
in Figs. H and |]. The soundness of the method 
of estimation of the coupling constants J\ and J2 
is illustrated in Fig. || where the helicity modulus 
is represented by the dashed line. Following the 
suggestion that J2/J1 can be estimated from the 
ratio between the dotted vertical lines we again 
get J2/J1 ~ 0.5. The quantitative differences be- 
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Figure 2. Re[wcr(a;)] vs T measured for 
a thin YBCO film at the same frequencies 
as Im[— uxr(u>)] in Fig.l. The inset illus- 
trates the peak ratio defined as the ratio 
Re[wcr(w)]/Im[— luo{lu)] taken at the dissipation 
peak maxima for a given lu. The data are taken 
from Ref. §. 



tween the data in Fig. g and the simulations in 
Fig. H is according to our interpretation due to the 
fact that the simulations cannot reach as low fre- 
quencies as the experiments (more precisely the 
ratio between the frequency and the microscopic 
frequency scale is larger in the simulations). A 
higher frequency broadens and smears the dissi- 
pation peaks. Since the interplane coupling in ex- 
periments presumably does not vanish, we have 
also investigated the influence of a small Jj_ in our 
model and found that a weak interplane coupling 
(0 < J± < 0.04) does not influence the qualitative 
features found above for J± = 0. 

The similarity with the experimental data by 
Festin et al. implies that the YBCO film is split 
parallel to the surface into two superconducting 
parts both having transitions with 2D character 
but with different T c . Since the coupling between 
the planes does not wash away the double peak 
transition one must infer that the coupling be- 
tween the two sheets is very weak. One possible 
reason for the phase separation may lie in a com- 
bination of an inhomogeneous surface with a more 
homogeneous part closer to the substrate. This 
can lead to different oxygen contents in the two 
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Figure 3. Im[— Lua(u>)] vs T for two weakly cou- 
pled 2D XY models with system size L = 64 at 
lu = 0,0.06,0.12,0.18 (from the left to the right). 
The coupling constants are chosen to be J% = 1, 
J2 = 0.5, and J± = 0. 



parts, which in turn results in slightly different 
lattice parameters causing a physical boundary 
between the two parts (|,[l^-[l|| 

A further consistency check on the scenario in 
terms of two coupled 2D superconducting parts 
comes from the peak ratio defined as the ratio 
Re[a;<7(w)]/Im[— wa(u>)] taken at the dissipation 
peak maxima for a given lu. For a 2D super- 
conductor this peak ratio should vary between 
2/n w 0.63 for small lu to 1 for larger lu Jl5[ . 
The inset in Fig. || gives the ratios estimated from 
the experimental data. The filled circles is for 
the large-T peaks and these ratios are consistent 
with 2D vortex fluctuations close to the transi- 
tion. The empty circles corresponds to the small- 
T peaks where in accordance with an interpreta- 
tion in terms of two transitions the zero level of 
Im[-w(r(w)] for the low-T transition is estimated 
by a linear extrapolation of the large T-part to- 
wards lower T. An almost as good estimate is to 
approximate by the plateau between the transi- 
tions (horizontal line in Fig. [l]) . These peak ratios 
are also rather consistent with a 2D transition in 
the limit of small lu (peak ratios for the lowest-T 
values). The small deviation for higher temper- 
atures may be due to a convolution of the high 
and low T parts of the transition. 
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Figure 4. Re[wcr(w)] as a function of temperature 
for different u): the lines with filled circles, empty 
circles and triangles correspond to the frequency 
0.06,0.12, and 0.18, respectively. Vertical lines 
show T c for two decoupled 2D XY models with 
coupling constants J\ = 1 and J2 = 0.5: T c \ w 
0.89 and T cl « 0.45. 



Another question is at what thickness d the di- 
vision takes place. A rough estimate may be ob- 
tained by assuming that both sheets consist of 
identical material. In such a case J\ oc d\ns and 
J2 oc d2ns where e?i(2) is the thickness of respec- 
tive sheets and n$ is the density of Cooper pairs 
for the material. Thus J2/J1 ~ cfe/^i and since 
we have found that J2/J1 ~ 0.5 from the data, 
we conclude that the boundary between the two 
parts should occur somewhere in the middle of 
the sample. 

In summary we conclude from our numerical 
simulations of two weakly coupled XY models 
that the double peak dissipation structure ob- 
served by Festin et al. is consistent with the 
interpretation that the sample consists of two 
parts with slightly different transition tempera- 
tures, and that these two parts are separated by 
a boundary which mainly runs parallel to the sub- 
strate. 
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